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PLASTIC BEHAVIOR OF CIRCULAR PLATES 
UNDER TRANSVERSE IMPULSE LOADINGS 
OF GAUSSIAN DISTRIBUTION* 

By Robert G. Thomson 
Langley Research Center 

SUMMARY 

An analytical study of the plastic response of a circular plate of uniform thickness 
subjected to a normal axisymmetric impulse of Gaussian distribution is presented. The 
Gaussian impulse loading is assumed initially to impart a momentum to the plate. The 
plate is thereafter prescribed as free from external loads. The kinetic energy of motion 
is assumed to be dissipated in plastic deformation. The radial and circumferential 
strain and strain-rate distributions of the plastically deforming plate are determined as 
a function of time and space from derived deflection expressions. Only the bending 
action of the plate is taken into account and the plate material is assumed to behave as a 
perfectly plastic -rigid continuum. 

The possibility of a fracture or separation of the plate material occurring during 
deformation is discussed. Maximum values of strain and strain rate are shown to occur 
in the circumferential components, and a mode of failure by radial fractures is indicated. 

INTRODUCTION 

The impingement and fracture of the main wall of a space vehicle by meteoroids 
are of major concern to the spacecraft designer. Protective devices, such as meteoroid 
bumpers, have been proposed to shatter the impinging meteoroid and to lessen its pene- 
trating potential. It has been shown both theoretically and experimentally (refs. 1 and 2) 
that the resulting cloud of meteoroid and bumper debris passes on to the main wall of the 
space vehicle as an impulsive loading having a general Gaussian distribution. A study 
of the dynamic response of a plate under this impulsive type loading could well classify 
the mode or manner of fracture at failure recently observed experimentally on double 
sheet models of a bumper and main wall configuration. (See refs. 3 and 4.) 

*The information presented herein was included in a thesis submitted in partial ful- 
fillment of the requirements for the degree of Doctor of Philosophy, Virginia Polytechnic 
Institute, Blacksburg, Va., 1966. 



A thorough analysis of the behavior of an impulsively loaded plate should include 
not only the plastic but also the elastic response of the plate material. Because of the 
complex time -dependent properties of the problem, however, most analyses consider 
only the rigid-plastic response of the plate and neglect its elastic behavior entirely. 

This rigid-plastic assumption for plate behavior is adopted in the present paper. The 
energy dissipation process then becomes one of pure plastic deformation in which a plate, 
once given an initial kinetic energy of motion, dissipates energy in plastic deformation 
and acquires a permanently deformed shape when the motion terminates. Energy trans- 
formation due to rotary inertia and dissipation due to transverse shear displacements 
are neglected. 

The existing analyses of the dynamic rigid-plastic behavior of plates consider the 
loading function as constant over the loaded area (refs. 5 to 9). Although this assumption 
on loading simplifies the mathematical development of the problem, it greatly restricts 
its application. The present analysis, therefore, considers a more general representa- 
tion of the loading function. The loading is still considered to be axisymmetric with 
respect to the center of the plate, but its radial distribution is given a general Gaussian 
shape. 


SYMBOLS 

b radius of finite circular plate 

C constant of integration (see eq. (17)) 

2h plate thickness 

k r radial curvature rate 

k^ circumferential curvature rate 

M Q yield-moment resultant, <7 0 h^; maximum bending moment plate can sustain 

M r radial bending-moment resultant 

M n circumferential bending-moment resultant 

u 
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q load distribution 

Q r transverse shear stress resultant 

r,0,z radial, circumferential, and axial plate coordinates, respectively (see fig. 2) 

S standard deviation 

t time 

t* time of cessation of all plate motion after impact 

tj time at which hinge -circle radius has decreased to zero 

V = w(r,t) 

V 0 = w(0,0) 

w plate deflection in transverse direction (z-direction) 

w* permanent plate deflection for t = t* 

e r ,eQ radial and circumferential strain, respectively 

6(t) Dirac 6-function, f 6(t) dt = 1 

V _oo 

ju mass per unit area of plate material 

p(t) radius of hinge circle 

p(0) initial position of hinge circle at t = 0 

Ujj stress components 

a r ,a 0 ,a z principal radial, circumferential, and transverse stress, respectively 

a Q yield stress in simple tension or compression (assumed to be identical in 

magnitude) 
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function of time (see eq. (32)) 


<Mt) 


A dot over a symbol indicates partial differentiation with respect to time t. 

FUNDAMENTAL CONCEPTS 
General Considerations 

The present analysis determines the rigid-plastic bending deformations, strains, 
and strain rates experienced by a circular plate when subjected to an axisymmetric 
Gaussian impulse loading. The circular plate, of radius b, is taken to be simply sup- 
ported, uniform in thickness, and composed of a plate material that is incompressible, 
ideally plastic (no work-hardening), and isotropic in yielding. Only the bending action of 
the plate is taken into account and the elastic deformations as well as the membrane 
forces are neglected. The plate material behavior is therefore considered to be per- 
fectly plastic -rigid and is assumed to deform according to the Tresca yield condition and 
the associated flow rule. 

The structural loading on the plate is assumed to result from a momentum exchange 
with an axisymmetric momentum impulse having a general radial Gaussian distribution 

where S is the standard deviation of the Gaussian momen- 
tum distribution. (See fig. 1.) This momentum exchange 
can be considered as an impulse loading such that at t = 0, 
a Gaussian distribution of load is instantaneously imparted 
to the entire plate. Thereafter, the plate is considered to 
be free from external loads. 

Equation of Motion 

Let r, 0, and z be a cylindrical coordinate system whose origin is located at 
the center of the undeformed middle surface of the circular plate and whose Z-axis is 
pointing vertically downward. (See fig. 2.) It follows from the symmetry of loading that 
the only nonvanishing stresses are o r and a Q and the transverse shear force Q r . 
Integration of a r and Oq through the thickness of the plate leads to the moments M r 
and Mg. 

The forces and moments acting on an element of the plate middle surface are then 
as shown in figure 2. Positive directions of the bending moments and shearing forces 
are shown, and all stress resultants are specified per unit length of line element of the 
middle surface of the plate. Summation of the vertical forces and summation of the 
moments acting on the plate result in the following equation of motion: 
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The stress profile through the thickness of a plate in rigid-plastic theory is differ- 
ent from that postulated in the bending of an elastic plate. If the rigid -plastic plate is to 
bend, each horizontal layer through the 


plate thickness must be at the yield con- 
dition. This assumption, for the Tresca 
yield condition, creates a discontinuity 
of stress at the middle surface which is 
allowable within the framework of plastic - 
rigid theory (ref. 10). The resulting 
principal moments (employing Tresca 1 s 
yield condition) are thus written as 




z 


Figure 2.- Element of the circular plate with applied forces 
and moments. 


where 2h is the plate thickness, and a r and Og are positive constants for positive 
values of z. Hence, the yield moment M Q equals o Q h^, where a Q is the yield stress 
in simple tension or compression. 


Strain-Displacement Relations 

The customary linear strain-displacement relations for strain are retained and 
particles originally on a normal to the plate middle surface are considered to remain 
normal to the middle surface as the plate is bent. The rates of radial and circumfer- 
ential straining for the present axisymmetric problem can then be written as 
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(4) 


e-r = 


-z = Zk 

9r 2 Zk r 


e 9 ~ 


. z l9w, z k 

r dr 9 


(5) 


where k r and are the radial and circumferential curvature rates, respectively. 
Consequently, 

fr = K ( 

^ ~ k e 

Plasticity Relations 

The idealized perfectly plastic-rigid response of the plate material is character- 
ized by a constitutive equation. The dependence between the plastic strain-rate com- 
ponents e and the stress components is formulated according to the flow (or 

incremental) theory of plasticity. (See ref. 11.) For the case of ideal plasticity, the 

. p 

strain-rate components e ^ are expressed simply as 

(i,i= 1,2.3) 


where 

f yield function which depends only on stress state 

G scalar function depending, in general, on spatial coordinates and time; G deter- 
mines magnitude of plastic strain rates 

For ideal plasticity the plastic strain-rate components e?j are directed normal 
(outward direction positive) to the yield surface, f being constant. The yield sur- 

face considered in this paper is that corresponding to the criteria of Tresca. It is con- 
cluded from the Tresca yield condition that yielding occurs when the maximum shear 
stress reaches a certain value or, mathematically, when 

a l " a 3 = a o 

where cr^, anc * 03 are the principal stresses and 0 ^ = 0 2 = 03 . Tresca* s yield 
condition for this problem is shown graphically in figure 3 in a a r ,a^ plane (a z = 0). 
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The plastic strain-rate components lie normal to the 
yield surface except at the corners of the hexagon 
where the normal is not defined but is bounded by 
the normals to the intersecting adjacent straight 
sides. The Tresca yield hexagon, therefore, not only 
defines the relationship between cr r and Oq but 
also characterizes the flow mechanisms (or curva- 
ture rates) of deformation. (See refs. 5, 10, and 12.) 

For a simply supported circular plate and the 
type of loading considered in this paper, the plate 
cannot reach the flow limit without first becoming 
plastic at the center; hence, through symmetry, Figure 3 . Tresca yie | d hex agon. 

M r = M q - M Q at r = 0. At the simply supported 

edge of the plate, the moment M r must vanish, or M r = 0 at r = b. The two plastic 
regimes of Tresca' s yield hexagon satisfying these boundary conditions and the equation 
of motion (eq. (1)) are regimes (point) A and (line) AB (fig. 3). Thus, the plate is initially 
divided into two different plastic regimes; an inner central plastic regime (called 
regime A) in which 




and a second outer annulus (called regime AB) in which 



where k r and kg are the radial and circumferential curvature rates defined in equa- 
tions (4) and (5). Note that in the outer annulus, the restriction of zero radial curvature 
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rate (eq. (8)) establishes, at most, a linear relationship in r for the velocity in that 
regime. For the velocities considered, the entire plate is in a plastic condition when 
deformation commences. 


Deformation Modes 


At the interface between the plastic regimes A and AB, a hinge circle is formed 
(as indicated in fig. 4) because of a discontinuity in the rate of circumferential curva- 
ture between the two regimes. (See ref. 5.) It will be found subsequently that this 

hinge circle will have a maximum radius 
at t = 0 when the kinetic energy 
imparted to the plate is a maximum, and 
will decrease thereafter and move 
toward the center of the plate. When 
the hinge circle finally shrinks to zero, 
the plate motion does not cease, but the 
plate continues to deform in such a man- 
ner that an inverted cone is formed with 
the hinge at r = 0, until the velocity 
becomes zero everywhere. 

Certain higher order derivatives, 
when passing through the hinge circle, 
are found to be discontinuous according 
to the plasticity theory developed for the 
dynamics of plastic circular plates. 

(For details, see refs. 5, 8, and 13.) Since the deflection w(r,t) obviously must remain 
continuous across a hinge circle, 




Figure 4.- Mode of deformation. 


dw 

dt 


<>-[*] ♦*!!?] 


(9) 


where the brackets denote possible discontinuities in the enclosed quantities in passing 
through the hinge circle and p is the rate of change of the hinge circle with time. The 
velocity of the plate particles w and the slope of the plate 9w/9r must remain continu- 
ous for moving hinge circles. The following continuity expressions can also be shown to 
be continuous across r = p(t) even though they are composed of members (that is, 
Q^w/W^, w, and 9w/9r) which may be discontinuous: 


3w 


9 2 w 

_9r_ 

+ p 

_ar 2 _ 


( 10 ) 
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[w] ♦ p|f] * 0 


( 11 ) 


Although equations (10) and (11) are not used directly to obtain a solution for 
w(r,t), they are employed as a check on the validity of the radial and circumferential 
strain and strain rates derived from w(r,t). 

MATHEMATICAL ANALYSIS 

The analytical derivation of the deflection of the plastically deforming plate is 
divided into two sections. The first section presents the solution for the motion of the 
plate when the plate is divided into two distinct plasticity regimes, separated by a hinge 
circle, up to the time of the disappearance of this hinge circle. The second section con- 
tains the solution for the motion of the plate as it continues to deform, up to the time of 
cessation of all motion. From the deflection expressions, the velocities, strains, and 
strain rates can be determined and are evaluated in the appendix. 

Solution With Hinge Circle of Finite Radius 

Velocity determination.- For a perfectly plastic -rigid circular plate initially sub- 
jected to a general Gaussian distribution of momentum, the resulting velocity distribu- 
tions are as follows: 

In regime A, the expression for velocity is obtained by direct integration of equa- 
tion (1) with respect to time, it being noted that M r and M^ are both equal to M Q 
(eq. (7)), 

w = V Q e a r (0 ^ r s p(t)) (12) 

Note also from equation (12) that the curvature rates k r and k^, in regime A, are non- 
negative (eq. (7)). 

In regime AB, however, the determination of the proper velocity distribution is not 
as straightforward. The first condition that can be utilized is that the radial curvature 
rate k r must be zero (eq. (8)) and thus, the velocity has at most a linear dependence 
in r. The velocity must be continuous across the hinge circle, and from the boundary 
condition of zero displacement at the simply supported edge (eq. (2)), must have zero 
velocity at the simply supported edge. The satisfaction of these conditions leads to 
the following expression for the velocity distribution in regime AB: 
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(p(t) srsb) (13) 


w= V 0 e 


-a 2 [p(t)]‘ 


b 
b - 


- r 


p(o 


where p(t) is the radius of the hinge circle which is, as yet, undetermined. 

The equation of motion (eq. (1)) is therefore not uniquely satisfied in regime AB 
until p(t) is determined as a function of t. In order to solve for the proper time vari- 
ation of p(t), the equation of motion is utilized, together with the velocity expression 
(eq. (13)) and the associated boundary conditions for regime AB 


M r (p,t) = M 0 
M r (b,t) = 0 


(14) 


where 0 S M r ^ M 0 for p ^ r % b. 

Equation (13) is differentiated to obtain the following expression for the 
acceleration: 


w= V n e' a V -bjL* (l - 2a2 b p + 2 a 2 p 2 )p 
° (b - p) 2 V ' 


(15) 


which is substituted into equation (1) with Mg set equal to M Q (t > 0). Integration of 
equation (1) then determines an expression in M r (r,t). Applying the boundary conditions 
on M r (eq. (14)) determines the constant of integration and leads to the following first- 
order total differential equation in p(t): 


12M q 

pV 0 b 2 


dt = 



2a 2 bp + 


2a2p2^ 


-1 



( 16 ) 


Equation (16), when integrated, yields (from ref. 14) 


12M n 

^t+C 

pV 0 b 2 



-J_fi - e.) 

(ab) 2 l by 




3/? 

(ab) 3 


erf(ap) 


(17) 


The constant of integration C can be evaluated in terms of p(0), the initial loca- 
tion of the hinge circle. Thus for t = 0 and p(t) = p(0), C becomes 


C = 





2 p«>r 


erf(ap(0)) 

(ab) 3 


(18) 
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The initial location of the hinge circle radius p(0) is determined from two 
requirements: (1) that the time rates of curvature k r ,k 0 be nonnegative in regime A 

and (2) that dp/dt (eq. (16)) be negative (since the kinetic energy imparted to the plate 
is a maximum at t = 0 and decreases thereafter). Also, of course, the value deter- 
mined for p(0) must be unique. The value of p(0) satisfying these requirements is 
found to depend on the shape of the initial Gaussian momentum distribution. For a 
Gaussian distribution that resembles a localized load (ab > /sT), it is found that the 
requirement that dp/dt be negative governs the determination of p(0). Hence, solving 
equation (16) for the maximum value of p(0)/b that will still yield negative values of 
dp/dt results in 


£(21 = 

b 



(2 ^ ^2ab < °°) (19) 


For a Gaussian distribution that resembles a more distributed load (ab < fe), the 
governing criteria for the determination of p( 0) is that the rate of curvature be non- 
negative in regime A and also that the value of p(0) be unique. Thus, from equa- 
tions (4) and (12) the rate of radial curvature k r can be written as 


fc = = 2V 0 a 2 ( 1 - 2a2 r 2) e ' a2r2 (20) 

9r 2 

and for k r to remain nonnegative at r = p(0) and for p(0) to be unique for all values 
of ab, p(0)/b must be equal to 

fi(2l - 2 = 1 _ (l i /2ab i 2) (21) 

b b /2ab 

The lower limit on ab (ab = I/V2) places the hinge circle at the outer edge of the plate. 

For values of ab < l//2, the hinge circle remains at the outer edge of the plate 
since the boundary condition M r (b,t) = 0 (from eq. (2)) prohibits regime A from encom- 
passing the edge of the plate and passing beyond. Hence, for ab < l/j2, 

£(21 =1 (0 < \f2ab § l) (22) 

b 

The lower limit on ab may only approach zero since the integral evaluation of 
equation (17) is invalid when a = 0. However, solutions do exist in the neighborhood of 
zero and are sufficient for practical application. The lower limit on ab (ab - 0) repre- 
sents a uniformly distributed impulse over the entire plate. 
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The value of t at the time the hinge circle vanishes, or t = t^, can be deter- 


mined from equation (17) by setting £ = 0 to yield 


12M r 


MV 0 b 


ti = -1 + 


2 


(ab)‘ 


- C = -1 + 


(ab) 2 


1 - 


«_(i . e!9l\ + fiM . 6 (effll) 


(ab)' 


\ b ) 




e -a p(0) _ er f(ap(0)) 
(ab)^ 


( 23 ) 


12M, 


The expression for tj is a cubic in ab and is plotted as a function of — ti in 

V 0 b 2 p 

figure 5. As is shown in figure 5, different expressions for p(0)/b are used in the 



Figure 5.- Time interval tj for initial hinge-circle radius to decrease to zero 
as a function of initial Gaussian distribution parameter ab. 


three different ranges of ab. (See eqs. (19), (21), and (22).) Note that as ab 

u 12M o x . 12M 

approaches °°, — t^ approaches 0 and as ab approaches 0, 2_ti 

MV 0 b^ p V 0 b 2 

, , ,, . , 12M n 

approaches 1; this value corresponds to the value for — tx presented in refer- 

MV 0 b 2 

ence 6 for a uniformly distributed impulse. 

Deflection .- The deflection of the plate can now be determined since the velocities 
of deformation have been established. The plate, during deformation, is divided into 
three distinct regions of interest as shown in figure 6. 
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For elements of the plate located between the 
outer edge of the plate and the initial position of the 
hinge-circle radius (p(0)), or in region I of figure 6, the 
deflection expression is dependent only on the velocity 
as given for the plastic regime AB (eq. (13)) and can 
be written as 


nr 

- pm 


n 


■ P (0) - 




Figure 6.- Hinge circle and plasticity 
regimes. 


w(r,t) = V ( 


£ 


,-a 2 P 2 


b - p 


dt 


(p( 0) ^ r £ b) (24) 


In region II of figure 6, elements of the plastic plate are first located in the plastic 
regime A and then as the hinge -circle radius passes over, in the plastic regime AB. The 
deflection is therefore dependent not only on the motion due to the velocity of the 
regime A when the point in question was in regime A, but also after the time t(r) (when 
the hinge circle passes through that point of p = r) on the deflection due to the velocity 
present in regime AB at that point up to the time t. 

w(r,t) = V 0 e _a2r2 t(r) + V 0 (p(t) £ r £ p(0)) (25) 

In region III, the elements of the plate are always situated in the plastic regime A. 
The deflection is therefore dependent only on the velocity present in regime A (eq. (12)) 
and can be written 


w(r,t) = V 0 e" a2r2 t (0 £ r £ p(t)) (26) 

Expressing dt in terms of dp (eq. (17)) and integrating equations (24) and (25) yield 
(from ref. 14) 


pV n b 

w(r,t) = V n — — 
’ ° 12M„ 


M) 


3/e\ 2 

1 2 \b / 


3 _ 1 _ £ 

2(ab) 2 2 b 


>-2a 2 p 2 


§W 




2(ab) 


2 2 


; -2a 2 p(0)2| 


(p(0) = r = b) (27) 
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where t(r) is found from equation (17). 

For t(r), p/b in equation (17) is replaced by r/b, because t(r) is the time it 
takes for the hinge circle to pass through the point p = r. Thus t(r), from equa- 
tions (17), (18), and (23) can be written as 



The time the hinge circle vanishes t^ is given by equation (23) and p(0)/b is defined 
over different ranges of ab by equations (19), (21), and (22). 

Solution With Hinge Circle of Zero Radius 

When the hinge -circle radius shrinks to zero at t = tj and the entire plate is in 
the plastic regime AB, the previous analysis must be terminated because the accelera- 
tion vanishes and the equation of motion reduces to (see (eq. (1)) 

J^(rM r ) - m o (30) 

or 

c 3 

M r = M o + — (31) 

Since M r = M 0 at r = 0, Cg is zero and M r remains at the constant yield moment 
value of M 0 for all radial values; this result is in contradiction to the prescribed plas- 
ticity condition in regime AB (see eq. (8)) and the simply supported edge condition 
(eq. (2)). 

Since the entire plate is now in the plastic regime AB, except for the origin, the 
radial curvature rate (eq. (4)) is zero. Integration of the radial-curvature -rate expres- 
sion, while the boundary condition of zero deflection at r = b and the continuity 
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conditions on the deflection and velocity at t = tj are maintained, yields the following 
expression for the plate deflection: 

w(r,t) = <p( t) ^1 - g) + w(r,t x ) (t g t x ) (32) 

where w(r,tj) is the deflection at t = tj, given by equations (27) and (28) with p set 
equal to zero, and where <p(t) is an arbitrary function. The plasticity condition of zero 

radial -curvature rate is satisfied if > o, This condition on cp( t) is examined for 

its validity once (p( t) has been determined. Note, however, that the slope of the deflec- 
tion as given by equation (32) is not zero at r = 0, because of the zero radial -curvature- 
rate condition. 

By a procedure similar to that used in the previous section of the paper, the accel- 
eration is obtained by differentiation of equation (32) and is substituted into equation (1) 
with Mq = M 0 . The subsequent integration of equation (1) and application of the boundary 
conditions 


yield 


M r = M 0 (r = 0) 

M r = 0 (r = b) 



(33) 


(34) 


which upon integration, by utilizing the conditions evaluated at t = tj 

<H n) = 0 1 

> (35) 

^(tl)= V o J 

becomes 

<#> ( *) - - - ‘l) 2 + V 0 (t - t t ) (36) 

It can be seen from equation (36) that 4> is indeed greater than or equal to zero, and 
w(r,t) for t > tj can be written as 
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(37) 


w (r,t) = w(r,tj) + (l - g) 


6M 


jtib 


a 


2V 'll +V o(< 


t - t 1 



The time t = t* at which the entire plate comes to rest can be found by differentiating 
equation (37) to obtain the velocity, by equating this expression for velocity to zero, and 
by solving for t*. Hence, 


12M 0 

MV 0 b2 


1 + 


12M 0 

fiV Q b 2 


(38) 


and the permanent deflection shape, denoted as w*(r,t*) for t = t* can then be written 
as 


w*,r,t*)»w ( r,t l)+ (l-£)^^ 


(39) 


The complete time -dependent expressions for the deflection shape of the plastically 
deforming plate, beginning with the initiation of motion, continuing past the time t = t^ 
(the vanishing of the hinge circle), and continuing up to the time of cessation of all plate 
motion t = t* have thus been determined. The complete time -dependent strains and 
strain rates of the plastically deforming plate can then be determined and are so derived 
in the appendix. 


RESULTS AND DISCUSSION 

In the evaluation of the permanently deformed shape of the plastic plate as the 
velocity drops to zero, it has been tacitly assumed that the plate remains continuous and 
homogeneous. In reality, however, "separation'* or "fracture" of the plate material could 
occur, under certain loading conditions, if the plate experiences "critical" strain or strain 
rates or combinations of these over portions of its surface. Under these conditions, of 
course, the application of the analysis should be terminated. After a choice of criteria 
is made, the knowledge of the complete time -dependent strain rate, strain, and velocity 
histories experienced by the Gaussian impacted plate permits the analyst to determine 
whether the integrity of the plate material has been maintained. In addition, a study of 
the strain and strain rate distributions can indicate the manner or mode of separation or 
fracture of the plate material when failure does occur. For these reasons, two sample 
Gaussian impulse loadings have been considered, and the velocity, strain rate, and strain 
distributions have been determined at various times after impact. Radial distributions 
of the velocities, strain rates, and strains are plotted for the two sample loadings in 
figures 7 to 11. 
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(d) ab = 2.0, ti < t < t 


for two sample loadings. 




12M 

_B 5- t 

plO)' k 2 




12M 0 12M 0 

(al ab = 1.0, 0 < t < ti; ^ ti = 0.467. (b) ab = 2.0, 0 < t < t 1; - = 0.074. 

MVgb^ MVob^ 

Figure 8.- History of radial strain rate for two sample loadings. Radial strain rate is zero after hinge circle vanishes, t > tj. 


In parts (a) and (b) of figures 7 to 11 the times chosen range between t = 0 and 
t = tj, during which time interval the hinge circle has a finite radius. The specific loca- 
tion of the hinge radius for the particular times chosen 0-8, 0-6, 0.4, and 0.2j 

is indicated in the figures by a circle. In parts (c) and (d) of figures 7, 9, 10, and 11 the 
times chosen range between t = t^ and t = t*. In this time interval, the hinge circle 
has already shrunk to zero radius; thus, the plate is entirely within the plastic regime AB. 
Parts (a) and (c) are calculated by using a value of ab of 1.0 and parts (b) and (d) are 
calculated by using a value of ab of 2.0. These sample values of ab (that is, Gaussian 
shape parameters) were chosen to fall within different ranges of ab for which different 
expressions are valid for p(0)/b. (See fig. 5.) 

Velocities 

The velocities are shown plotted in figure 7. The velocity as shown in parts (a) 
and (b) exhibits Gaussian character in the central region of the plate and decreases 
linearly with increasing radius in the surrounding annulus. The central region is seen 
to decrease progressively as time increases (p decreasing) until the entire plate is 
within the plastic regime AB and the velocity is completely linear in character. The 
velocity as represented by parts (c) and (d) of figure 7 remains linear, the magnitudes 
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decreasing as the time increases. These velocity profiles are as expected, since the 
basic plasticity assumptions made with respect to the rates of curvature only admit of 
these profiles in the two plastic regimes (eqs. (12) and (13)). 

Radial Strain Rate 

The radial strain rates are plotted in figure 8. The radial strain rates are zero in 
plastic regime AB because of the basic plasticity assumption made with regard to the 
rates of curvature possible in this regime (see eqs. (4) and (8)) and indeed the radial 
strain rate (as given by eqs. (A6) and (A18)) is zero for all times. In figures 8(a) and 
8(b) the radial strain rate possesses an approximate Gaussian radial distribution in the 
central region of the plate (regime A) with a discontinuous decrease to zero across the 
hinge circle to the surrounding annulus of regime AB. Note that the strain rate does not 
vary with time but the movement of the hinge circle toward the center of the plate 
decreases the central plastic region (regime A) and causes a decrease in the percentage 
of the plastic plate with nonzero radial strain rate. The magnitude of the discontinuity 
present at the hinge circle, however, increases with decreasing hinge radius. 

Circumferential Strain Rate 

The circumferential strain rates are plotted in figure 9. In parts (a) and (b) of 
figure 9, the circumferential strain rates also exhibit a Gaussian radial distribution of 
the central region (plastic regime A) of the plastic plate; however, a large increase in 
magnitude occurs as one passes into the surrounding annulus of regime AB. The circum- 
ferential strain rate in regime AB is the maximum value of strain rate experienced by 
the plastically deforming plate. The circumferential strain rates in regime AB are 
inversely proportional to the radius and decrease accordingly as the radius increases. 

The magnitude of the discontinuity across the hinge circle is similarly seen to increase 
as the radius of the hinge circle decreases. In parts (c) and (d) of figure 9, the entire 
plate is contained within the plastic regime AB and the circumferential rates are shown 
to decrease in magnitude with increasing time with a radial distribution again propor- 
tional to the inverse of the radius. 


Radial Strain 

The radial strains are shown plotted in figure 10. The central region of the plastic 

plate exhibits an approximate Gaussian type radial distribution of radial strain (plastic 

regime A) with a discontinuous "jump" in strain developing across the moving hinge 

circle in passing into the plastic regime AB. This discontinuity in radial strain at the 

moving hinge circle is the jump condition noted in equation (10) for 9^w/9r^ (divided 

by -h) ; a second discontinuity develops at -= for ^ — (fig. 10(a)), but not 

3 b b b ^2ab 
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(fig. 10(b)). The plastic regime AB always has zero radial 
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strain. For regime A, however, as the hinge circle moves toward the center of the plate, 
this central plastic region decreases and the surrounding annulus of regime AB encom- 
passes regions of the plate formerly situated in regime A. This annulus thus contains 
residual strains from regime A, as well as the radial straining that occurs across the 
discontinuity of the moving hinge circle as is shown in parts (a) and (b) of figure 10. As 
is to be expected, the radial strains begin at zero and increase in magnitude as the time 
progresses. In parts (c) and (d) of figure 10, the residual radial strains present in the 
plastic regime AB (caused by radial straining in the plastic regime A and radial straining 
at the discontinuity of the hinge circle) remain constant with time, for t > t^, and no 
further increase in radial strain occurs because the radial straining in regime AB is 
zero. 


Circumferential Strain 

The circumferential strains are plotted in figure 11. In parts (a) and (b) of fig- 
ure 11, the circumferential strains again exhibit Gaussian radial distribution in the inner 
central region (plastic regime A), a discontinuity developing in the slope of the circum- 
ferential strain as one passes through the moving hinge circle into the surrounding 
annulus. In the surrounding annulus, the magnitudes of the circumferential strains 
increase since this region contains both residual strains developed when this part of the 
plate was in the plastic regime A and further circumferential straining when this part 
passes into the plastic regime AB. A second discontinuity in slope occurs when the 
initial position of the hinge circle is reached if the initial position is determined from 
P,(9) = _J_. (See fig. 11(a).) The circumferential strains begin at zero and increase in 
k t/2ab 

magnitude as the time increases. In parts (c) and (d) of figure 11, the residual circum- 
ferential strains developed in regimes A and AB are present (and remain constant with 
respect to time) and additional circumferential straining is experienced by the plate (now 
completely in regime AB) as time progresses. The circumferential strain continues to 
increase in magnitude as time progresses and reaches a maximum as the impacted plate 
comes to rest. (Note again that in the immediate neighborhood of r = 0, the solutions 
presented are invalid.) 

The increase in magnitude of the circumferential strains and strain rates as one 
passes from the central region (plastic regime A) of the plate to the surrounding annulus 
(plastic regime AB) could have a special significance. With the passing of time, the 
hinge circle decreases in radius, as has been noted before, and the circumferential 
strain rates and strains in the plasticity regime AB immediately adjacent to the hinge 
circle subsequently increase in magnitude. This increase in the magnitude of the 
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circumferential strain rate and strain eventually exceeds that of their radial counterparts 
as the hinge -circle radius shrinks to zero. If failure by fracture or separation does 
occur during the deformation process, the magnitudes of the circumferential components 
probably would be the maximum values of strain rate and strain experienced by the plate 
at failure, and fracture would occur along radial paths. 

CONCLUDING REMARKS 

An exact analytical solution is presented for the plastic behavior of a circular 
plate subjected to a general radial Gaussian impulse loading. Axial symmetry being 
assumed, a small deflection bending analysis is developed to determine the complete 
time -dependent strain and strain-rate distributions present in the plastically deforming 
plate from the moment of impact (considered to be instantaneous) to the cessation of 
motion. Graphical representations of the radial and circumferential components of the 
strain and strain rate are presented for various times between initial impact and 
cessation. 

A study of these graphical representations indicates that the maximum values of 
strain rate occur in the circumferential direction and that the radial strains exhibit a 
discontinuity in crossing through the moving hinge circle which separates the deforming 
plate into two distinct plasticity regimes. The maximum values of strain and strain rate 
are also shown to be distributed initially in annular rings about the center of the plate. 
These maximum values increase in magnitude as time progresses, and the hinge circle 
moves toward the center of the plate. The circumferential components eventually 
increase to become the maximum values in both strain and strain rate experienced by 
the plate. 

A study of these strain and strain-rate distributions suggests a pattern of fracture 
of the plate material when failure occurs. Because the circumferential strains and 
strain rates are higher than their radial counterparts, fractures can be expected to occur 
along radial lines. These radial fractures may initiate the catastrophic pe tailing -type 
failure noted experimentally on the inner wall of double-sheet micrometeoroid protective 
devices. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 22, 1967, 

124-08-01-13-23. 


24 



APPENDIX 


VELOCITIES, STRAIN RATES, AND STRAINS OF 
THE PLASTICALLY DEFORMING PLATE 

The complete time -dependent expressions for the deflection shape of the plastically 
deforming plate have been determined in the body of the paper. In this appendix, the com- 
plete time -dependent velocities, strain rates, and strains are derived. 

Solutions for Times Prior to Shrinking 
of Hinge Circle to Zero Radius 

Velocities.- The velocities are as given in equations (12) and (13) or can be derived 
by differentiation of equations (24) to (26) with respect to t, t(r) being considered con- 
stant with respect to time, 


w(r,t) = V 0 e _a ^ r ^ 


(0 S r i p) (Al) 


w(r,t) = V Q e -a P 


2 n 2 b - r 


b - P 


(p ^ r S b) (A2) 


Strain rates.- The strain rates, as given by equations (4) and (5), are linear in z, 
the maximum strain rate occurring at z = ±h. Consider the maximum strain rate in 
tension for a positive bending moment (z = h) (from eqs. (4) and (5)): 


h 8 Vr,t) 
0 

(A3) 

dr* 


h 9w(r,t) 
r Or 

(A4) 


Differentiation of equations (Al) and (A2) with respect to r accordingly yields for the 
radial strain rate 


e r = 2hV 0 a2e“ a ^ r ^ 

CM 

CM 

1 

1 


(0 = r H p) 

(A5) 

e r = 0 


"D 

IIA 

»-s 

IIA 

cr 

(A6) 
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and for the circumferential strain rate 


= 2hV 0 a^e“ a ^ r ^ 


(0 g r ? p) (A7) 


^A Ve -aV 1 

6 br 0 




(p = r g b) (A8) 


Strains .- The maximum strain in tension for a positive bending moment (z = h) 
from equations (A3) and (A4) can be written as 


= _ h 9 w ( r >t) 


8r" 


(A9) 


_ h 8w(r,t) 
“ r 9r 


(A10) 


Differentiation of equations (24) to (26) with respect to r, accordingly, yields for the 
radial strain 


e r (r,t) = 2ha 2 V 0 e _a2r2 




(0 r 5p) (All) 


e r (r,t) = 2ha 2 V 0 e“ a2r2 


1 - 2< ab > 2 (£) 2 


t(r) + A ' ,V ° b e -2a 2 r 2 
b 2 12M 0 


1 - 2(ab) 2 g 


+ «-» ^) 2 


n 2 


(P S r S p(0)) (A12) 


e r (r,t) = 0 


(p(0) Ir ^b) (A13) 


Note that t is the time when the hinge circle reaches the radial position r and 
that t(r) is the time the hinge circle crosses the point p - r. (See eq. (29).) Simi- 
larly, differentiation of equations (24) to (26) with respect to r for the circumferential 
strain yields 


e 0 (r,t) = 2hV 0 a 2 e" a2r ^t 


(0 5r5p) (A14) 
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9 c. 2,.2 , , V n h pV n b z 

£ „(r,t) = 2hV 0 a^e -a t(r) + — — 

b 


3 /£\ 2 

2 V b / 2 (ab ) 2 2 b _ 


_3 1 _ £ 


e - 2 a 2 p 2 


-f- ) + - 

2 \b / 2 (ab) 


1 r 
2 " 2 " b 


. - 2 a 2 r 2 


(p = r ^ p(0)) (A15) 


C/j(r,t) = 


V 0 h pV Q b 

b 2 r 12M 0 
b 


sfeM 

2 \ b 


2 [ r; 


3/a 

2 (b 


3 _ 1 _ £ 

2 (ab ) 2 2 b 


e -2aV 


3 1 P(0) 


2 (ab) 


2 2b 


1 - 2 a 2 p( 0 ) 2 | (p(o) lr?b) 


(A16) 


Solutions for Times After Disappearance of Hinge Circle 

The velocities, strain rates, and strains for tj £ t ^ t* can be obtained by the 
use of equation (37). 

Velocities.- Differentiation of equation (37) with respect to t yields 


v(r,t) = V Q (b - r) 


w 


12M 


M V 0 b 


3(*» - ‘) 


+ 1 


(0 S r s b) (A17) 


Strain rates. - The maximum strain rates in tension as found from equations (A3) 


and (A4), by differentiation with respect to r of equation (37), become 


e r = 0 


(0 ? r 5 b) (A18) 


e 9 = 


hV 0 

b 2 I 


12M 


b •- 


,iV 0 b 


3(h • l ) + 1 


(0 < r ? b) (A19) 


Strains.- The maximum strains in tension for a positive bending moment (z = h) 
are found from equations (A9) and (A10) by the proper differentiation of equation (37) 
with respect to r, and can be written as follows: 
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e r (r,t) 




2e- a2 r 2 

1 - 2(ab ) 2 f-f 

t(r) | e' 2a2r2 V G pb 2 f 


\b j 

tl (ab) 2 12M 0 t l \ 


+ 2(ab) 2 ^ 


si 2 


1 + 3 


1 - 2(ab) 2 | 


(0 g r 5 p( 0)) (A20) 


e r ( r >t) 

hV 0 a 2 t! 


= 0 


(p(O)Sr^b) (A21) 


e lM = 2e -a 2 r 2 tfcl + 1 V o b 

hV 0 a 2 t 1 h (ab) 2 | 12 M Q t 1 ^ 2 (ab) 2 2 


_ 1+ 3/r ' 2 


2(ab) 2 2 2\b) b 


>-2a 2 r 2 ' 


+ 


(ab) 2 l 


6M o (> -tl) 2 ( t t 
pV Q b 2 tl H 


(0 < r S p(0)) (A22) 


e 0 (r,t) 


pV Q b^ 


f 3 1 

3/p(0)V 

Wab) 2 2 

2 \ b ) 


2(ab) 2 


,-2a 2 p(0) 2 ’ 


(ab) 2 | 


6M 


MV n b' 


o 1 . j 


(p(0) ^ r £ b) (A23) 


Note that in the neighborhood of -= 0, equations (A8), (A15), (A16), (A19), (A22), 

and (A23) are invalid. Since the principal rate of radial curvature is zero in regime AB 
(eq. (8)), a discontinuity in slope exists at r = 0. 

2 

The discontinuities in the second derivatives — ™ and w which appear as 


dr* 


one crosses from one side of the moving circle to the other can be shown to satisfy the 
discontinuity relationships (eqs. (10) and (11)). The velocity remains continuous, as was 
assumed initially (see eqs. (12) and (13)), and so does the radial slope 9w/9r (see 
eqs. (A10), (A14), and (A15)) in passing through the hinge circle. The deflection w, of 
course, also remains continuous. (See eqs. (26) and (28).) The discontinuity in 9 2 w/W 2 
is given by the second term in equation (A12) ^divided by -h and evaluated at g = 
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which is also the negative of the discontinuity in 9w/9r divided by dp/dt (eq. (16)). 

The discontinuity in 9w/9r is found from the difference between equations (A7) and (A8) 
when £■ divided by -h/r; this difference, according to equation (11), equals the neg- 
ative of the discontinuity in the acceleration divided by dp/dt. This result can be veri- 
fied with equations (15) and (16). 
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